Spectropolarimetry of the Peculiar Type la SN 2005hk 



Ryan Chornock 1 , Alexei V. Filippenko 1 , David Branch 2 , Ryan J. Foley 1 , Saurabh Jha 1 , and 

Weidong Li 1 



^ ABSTRACT 
O 

o . 

We present Keck spectropolarimetry of the unusual type la supernova (SN la) 
2005hk several days before maximum light. An analysis of the high signal-to- 
noise-ratio total-flux spectrum shows the object's extreme similarity to the pe- 
culiar SN 2002cx. SN 2005hk has an optical spectrum dominated by Fe III lines 
and only weak lines of intermediate-mass elements, unlike a normal SN la at this 
\ epoch. The photospheric velocity measured from the minima of strong absorp- 

tion lines is very low for an SN la (~6000 km s -1 ), solidifying the connection to 
SN 2002cx. The spectrum-synthesis code Synow was used to identify the pres- 
ence of iron-peak elements, intermediate-mass elements, and possibly unburned 



o 



carbon at similar velocities in the outer ejecta of SN 2005hk. Many weak spectral 
features remain unidentified. The spectropolarimetry shows a low level of con- 
tinuum polarization (~0.4%) after correction for the interstellar component and 
only a weak Fe III line feature is detected. The level of continuum polarization 
is normal for an SN la, implying that the unusual features of SN 2005hk cannot 
' be readily explained by large asymmetries. 
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Subject headings: supernovae: individual (SN 2005hk) — polarization 



1. Introduction 

Supernovae of type la (SNe la) represent the thermonuclear explosion of carbon-oxygen 
white dwarfs, probably as they approach the Chandrasekhar mass in a binary system. De- 
spite the importance of SNe la as sites of iron-peak nucleosynthesis and their utility as 
cosmological standard candles, the actual progenitor systems and details of the explosion 
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mechanism are still unknown; see Hillebrandt & Niemeyer (2000) for a recent review of 
theoretical models and problems associated with them. 

What is clear from an observational point of view is that most SNe la fall on a sequence 
between the peculiar SN 1991bg-like objects at one end and the peculiar SN 1991T-like 
objects at the other (Branch, Fisher, & Nugent 1993). The photometric sequence between 
these two endpoints (Phillips 1993) is paralleled by a spectroscopic sequence (Nugent et 
al. 1995). SN 1991bg (Filippenko et al. 1992b) was faint, red, had a rapidly declining light 
curve, and showed prominent lines of Ti II in its spectrum at maximum light while SN 1991T 
(Filippenko et al. 1992a; Phillips et al. 1992) was bright, blue, had a slowly declining light 
curve, and showed prominent Fe III lines in its spectrum. Nugent et al. (1995) explained 
the variation in normal SNe la between these two extrema as primarily a temperature effect, 
driven by the amount of radioactive 56 Ni formed in the explosion. Most objects are normal, 
but estimates of the peculiarity rate vary from 11% to 36% (Branch, Fisher, & Nugent 1993; 
Li et al. 2001a), depending on the definition of "peculiar." 

In recent years, some variants from this simple prescription for SN la diversity have been 
found (Li et al. 2001b, 2003; Garavini et al. 2005), including the intriguing object SN 2002cx. 
Li et al. (2003) showed that SN 2002cx had Fe III lines in its pre-maximum-light spectrum 
like SN 1991T, a low total luminosity comparable to SN 1991bg, the lowest photospheric 
velocities ever measured in an SN la at maximum light, an almost normal S-band light curve 
but a bizarre plateau-like light curve in R and /, peculiar color evolution, and unusually 
rapid spectral evolution. They concluded that no existing theoretical model could match 
the combination of features that SN 2002cx exhibited. Jha et al. (2006) showed that the 
peculiarities persisted to very late times. The spectrum of SN 2002cx nearly 300 days after 
explosion was dominated by a forest of narrow permitted Fe II lines, each with a velocity 
width of ~700 km s -1 . They also reported the detection of low- velocity intermediate- mass 
elements and possibly O I, with important consequences for SN la explosion models. In 
addition, Jha et al. (2006) presented spectra of several similar objects, including SN 2005hk, 
showing that SN 2002cx was not unique and therefore a new subclass of SNe la had been 
overlooked until recently. 

In this paper, we present Keck spectropolarimetric data of SN 2005hk, a peculiar SN la 
very similar to SN 2002cx, and an analysis of the total-flux spectrum. Full light curves and 
a spectral sequence will be presented elsewhere (Phillips et al. 2006, in prep.). 
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2. Observations 

SN 2005hk in UGC 272 was discovered (Burket & Li 2005) at an unfiltered magnitude 
of ~17.5 on 2005 October 30.25 (UT dates are used throughout this paper) by the Lick 
Observatory Supernova Search with the 0.76-m Katzman Automatic Imaging Telescope (Li 
et al. 2000; Filippenko et al. 2001; Filippenko 2005). An independent discovery on 2005 
October 28 was reported by the Sloan Digital Sky Survey II team, who measured a g mag- 
nitude of 18.8, two days after a non-detection (no upper limit was reported; Barentine et 
al. 2005). The host galaxy was observed as part of the Sloan Digital Sky Survey (SDSS 
J002749. 73-011200.0; Adelman-McCarthy et al. 2006) and a redshift of z = 0.0131 ± 0.0002 
was measured. We have adopted this redshift and removed it in all plots in this paper. 

Serduke et al. (2005) obtained an optical spectrum of SN 2005hk at Lick Observatory 
on November 2, which showed the object to be an SN la with prominent absorption lines of 
Fe III, similar to the peculiar object SN 1991T. The low photospheric velocities and faint 
apparent magnitude led Jha et al. (2006) to classify SN 2005hk as a member of an emerging 
subclass of SN 2002cx-like SNe la. 

We observed SN 2005hk on 2005 November 5.4 using the Low Resolution Imaging Spec- 
trometer (LRIS; Oke et al. 1995) on the Keck I 10-m telescope in spectropolarimetry mode 
(LRISp) 1 . We obtained a total exposure time of 5677 s in the four waveplate positions 
(0 = 0°, 45°, 22.5°, 67.5°). The unusual total exposure time was due to poor guiding during 
the first exposure of the set, which required us to abort and restart the waveplate rotation 
sequence. After the end of the sequence, a short (200 s) exposure in the second waveplate 
position (9 = 45°) was taken to complement the aborted first exposure. These two short- 
ened exposures did not substantially change the polarimetry if they were included in the 
reductions, so we have used them in the analysis below. 

Our instrumental setup employed both the blue and red sides of LRISp, with the D560 
dichroic splitting the incoming light beam. The combination of the 600/5000 grism on the 
blue side and the 400/8500 grating on the red allowed us to cover a wavelength range of 
3170 - 9240 A. The 1.5"- wide slit gave a spectral resolution of 6 A and 9 A on the blue and 
red halves of the spectrum, respectively. Conditions were good and the seeing was about 
0.9" during the SN 2005hk observations. 

Basic two-dimensional image processing was accomplished using IRAF 2 , along with 



1 See http://www2.keck.hawaii.edu/inst/lris/manuals/polarization/pol_v3.ps for the online polarimeter 
manual (Cohen 1996). 

2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
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optimal extraction and wavelength calibration of the one-dimensional spectra. Our own 
IDL procedures (Matheson et al. 2000) were used to flux-calibrate the spectra and correct 
for telluric absorption. The starting airmass was low (1.07), but it climbed to 1.2 by the 
end of observations. The slit was oriented at a position angle (P.A.) of 130° to minimize 
contamination from the host galaxy in background regions along the slit. The parallactic 
angle over the course of the observations varied from 10° to 50°, so differential light loss 
due to atmospheric dispersion (Filippenko 1982) is visible in the individual spectra of the 
polarimetric sequence. We assumed the first complete observation had the correct spectral 
shape, as it was taken at the lowest airmass and the P.A. was closest to parallactic. The 
spectral shape of the combined flux spectrum was corrected to that of the first observation 
by means of a low-order spline fit to the ratio of the two. 

The spectropolarimetric reductions followed the basic procedure of Miller et al. (1988), 
as implemented by Leonard et al. (2001). The polarimetric standard star HD 204827 
(Schmidt et al. 1992) was observed to determine the P.A. offset of the instrument from 
the sky reference frame. The blue (red) side flux-weighted average polarization angle over 
the wavelength range 3980-4920 A (5890-7270 A) was set to the catalogued value of 58.2° 
(59.1°) for the B (R) band. The derived offsets agreed to within 0.5° with observations of the 
polarization standard HD 19820 (Schmidt et al. 1992). Two null standards, BD +32 3739 
and HD 57702 (Schmidt et al. 1992; Mathewson & Ford 1970), were also observed and showed 
polarizations of <0.08%. 



3. Results 

3.1. Total-Flux Spectrum 

A by-product of the long exposure times necessary for spectropolarimetry is the high 
signal-to- noise-ratio total-flux spectrum. In Figure 1, the spectrum of SN 2005hk is compared 
to early-time optical spectra of several other SNe la. The SN 2002cx spectrum is from four 
days before maximum light in the B band (Li et al. 2003), about the same as this observation 
of SN 2005hk (W. Li private com.; Phillips et al. 2006, in prep.), and appears very similar. 

The most prominent features in the spectra of both objects are absorptions due to Fe III 
AA4404, 5129. These Fe III lines are one of the distinguishing characteristics of the peculiar 
SN 1991T-like class of objects, which are represented in Figure 1 by both the prototype, 
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SN 1991T (Filippenko et al. 1992a), and the well-observed SN 1997br (Li et al. 1999). One 
distinction is that the Fe III absorption minima in SN 2002cx and SN 2005hk are observed 
to fall near 4320 A and 5030 A, implying unusually low expansion velocities (~6000 km s _1 ) 
for these SNe la near maximum light (Li et al. 2003; Branch et al. 2004). 

Normal SNe la at this epoch also have prominent spectral features due to intermediate- 
mass elements (Si II, Ca II, S II, etc.), as marked in Figure 1 next to the spectrum of 
SN 1994D (Filippenko 1997). These features are much weaker in the other objects, but are 
clearly stronger in SN 2005hk than in SN 1991T. They are also stronger in this spectrum of 
SN 2005hk than in the Lick spectrum from three days earlier presented by Jha et al. (2006), 
indicating rapid spectral evolution, as was observed in SN 2002cx at early times (Li et al. 
2003). 

We used the parameterized supernova spectrum-synthesis code Synow (Fisher et al. 
1997; Fisher 1999) to investigate the SN 2005hk spectrum more quantitatively. Synow 
uses several simplifying assumptions for computational ease, such as spherical symmetry, 
a blackbody-emitting photosphere, and a resonance-scattering line source function in the 
Sobolev approximation. Each ionic species is assumed to have an optical depth that declines 
exponentially with velocity above some photospheric velocity, v p h, with an e-folding scale 
v e , and level populations in local thermodynamic equilibrium (LTE) at some excitation 
temperature (T exc ). The strength of Synow is that of a fast line-identification tool, rather 
than ejecta abundance determinations or detailed spectral modeling of the full SN spectrum. 

In the top panel of Figure 2, we present a fit using seven ions that we believe are present 
in SN 2005hk with high confidence. The photospheric velocity used in the fit was 5000 km s _1 
and the ions were assumed to have v e = 2000 km s -1 , with no maximum velocity imposed. 
(However, O I and Mg II are mildly detached, v p h = 6000 km s _1 .) The low-ionization species 
(Ca II, Mg II, S II, Si II, and O I) were fit using an excitation temperature of 10,000 K, while 
a higher T exc of 15,000 K was used for the doubly ionized species (Fe III and Si III). The 
higher value of T exc was necessary to match the strength of the Fe III feature near 5900 A. 
This synthetic spectrum fails to fit some of the features of SN 2005hk, particularly at bluer 
wavelengths, so we generated a second fit with four additional ions (Ti II, Ni II, Co II, 
and C III), whose presence we regard as possible, but not definite. This second synthetic 
spectrum is shown in the bottom panel of Figure 2. The additional ions are responsible 
for spectral features at desirable wavelengths, as well as a decrease in the amount of near- 
ultraviolet (UV) flux, as discussed below. A list of reference lines used in the fit and their 
optical depth can be found in Table 1. 

Branch et al. (2004) presented a Synow analysis of the SN 2002cx spectra. Of most 
relevance to us as a basis for comparison is the earliest spectrum they analyzed, from one 
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day before maximum light. The best fit for SN 2002cx has v p h = 7000 km s -1 and v e = 1000 
km s" 1 . While the derived parameters imply a somewhat higher photospheric velocity in 
SN 2002cx, Vph and v e are covariant and the two spectra show very similar absorption minima, 
except that in direct comparison it is evident that the minimum of the Si II A6355 line is at 
a lower velocity in SN 2005hk. The SN 2002cx Synow fit used six species (Si II, Si III, S II, 
Ca II, Ti II, and Fe III), all of which are present in our SN 2005hk fit. The presence of Ti II 
is not certain in either fit, but Ti II A3760 matches the absorption seen near 3680 A in both 
objects. Branch et al. (2004) observed that this SN 2002cx fit could provide an adequate 
model for the SN 1991T-like object SN 1997br (Li et al. 1999) if the photospheric velocity 
were increased, indicating the kinship between these two classes of supernovae. 

The SN 2005hk fit has five species not used for SN 2002cx. Two of the ions, O I and 
Mg II, have their strongest effects on the spectrum between 7500 A and 9500 A, so they 
could not have been detected at early times in SN 2002cx because the available spectra do 
not extend as far to the infrared (IR). The Mg II A4481 line also makes a contribution in the 
blue, but it could have been overlooked in SN 2002cx. We also note that the Mg II and O I 
lines in the red part of the spectrum are stronger than the neighboring Ca IR triplet, which 
is opposite of most SNe la (Filippenko 1997). By comparison, the early spectra of SN 1991T 
are rather featureless at these wavelengths. Oxygen and magnesium are both products of 
incomplete carbon burning, so their presence could be a diagnostic of the explosion process. 
We cannot be certain that the difference here between SN 2005hk and SN 1991T is solely 
due to differing composition as opposed to greater velocity smearing or higher temperatures 
and ionization levels, but more detailed study is warranted. 

In addition, Ni II and Co II are used in the second Synow fit shown in Figure 2, 
but were not used in the SN 2002cx fit. Their presence helps to fit features in the UV 
that were not included in the spectral range of the early SN 2002cx spectra. In particular, 
the flux minimum near 3300 A is rather interesting. SN 1991T showed a deep absorption 
trough near 3200 A that was interpreted by Filippenko et al. (1992a) to be the signature of 
Co II, and hence of complete silicon burning in the outer envelope of that object. Jeffery 
et al. (1992) showed that while iron-peak elements are responsible for the UV opacity in 
SNe la, an atmosphere consisting solely of iron would have a flux peak here, while nickel and 
cobalt produce a minimum. The more detailed models of Mazzali et al. (1995) confirmed 
and extended this result. By varying the composition of the atmosphere, they found that 
deflagration products (taken from the popular W7 model of Nomoto et al. (1984)) were 
inadequate to fit the depth of this feature in SN 199 IT at early times and a substantial 
amount of 56 Ni and its decay products was necessary. Synow is the wrong tool to use 
to determine abundances in SN 2005hk, but the identification of Ni II and Co II in this 
spectrum is suggestive of the presence of completely burned material in the outer layers of 
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SN 2005hk. 

The presence or absence of carbon in the spectrum is particularly interesting because 
of the implications for modeling the explosion of SN 2005hk. Some three-dimensional simu- 
lations of SN la deflagrations have shown that a significant amount (~40%) of the carbon 
and oxygen remains after explosion (Travaglio et al. 2004). One potential explanation for 
both the low luminosity and low expansion velocities of SNe 2002cx and 2005hk is that the 
explosions were very inefficient and an anomalously large amount of fuel was left unburned. 
At the conditions expected in an early-time SN la atmosphere, C I through C III are the 
relevant ionization stages of carbon (Hatano et al. 1999). Features attributed to C II A6580 
and C II A7235 have been identified in the pre-maximum spectra of several normal SNe la 
(e.g., SN 1998aq, Branch et al. 2003) and C I has been identified in the near-IR spectra 
of the peculiar low-luminosity SN 1999by (Hoflich et al. 2002), though Marion et al. (2006) 
placed strict upper limits on the presence of C I in the pre-maximum near-IR spectra of 
three normal SNe la. Unfortunately, our spectra do not extend far enough to the near-IR to 
test for the presence of C I lines. 

In the Synow fit presented in the bottom panel of Figure 2, we have included C III 
because an absorption near 4564 A is consistent with the strongest optical line of C III, the 
triplet at A4647, blueshifted by the photospheric velocity. A similar feature was tentatively 
identified as C III in an early spectrum of SN 1999aa by Garavini et al. (2004). We regard 
this detection as possible because we cannot claim an identification based on a single spectral 
feature to be unambiguous, especially given the many as-yet unidentified lines. The Synow 
fit shown in Figure 2 does show an additional C III feature just outside of our wavelength 
range at 9300 A, which might be observable in other spectra or other SN 2002cx-like objects. 

SN 2005hk has a weak emission feature at 7250 A with a corresponding P-Cygni ab- 
sorption at 7100 A that seems promising for C II A7235, but attempts to fit it using Synow 
failed. The Synow fits to A7235 overpredict absorption due to A6580. The observed spec- 
trum does have an absorption feature at the appropriate wavelength, but it is much too 
weak to match the model. LTE models generically predict that A6580 should be stronger 
than A7235 because the lower levels of A7235 are the upper levels for A6580. Likewise, the 
synthetic spectra with C II predict an absorption due to the A4745 multiplet which falls on 
an observed emission feature at 4662 A. The two longer-wavelength features are either not 
due to C II or they are far out of LTE, but in any case we cannot identify C II with certainty. 
Jeffery et al. (1992) encountered similar difficulties with the identification of C II in the early 
spectra of the more normal SN la 1990N (but see Mazzali 2001). The wavelength coincidence 
of features in the observed spectrum with the strongest expected lines of C II and C III is 
suggestive, but we cannot claim a firm detection of any form of carbon at this time. 
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One last unusual aspect of the SN 2005hk spectrum is the presence of many low- 
amplitude, high-frequency spectral features, particularly at wavelengths less than ~6000 A. 
See Figure 3 for a closer look at some of these features, which have amplitudes of a few 
percent of the continuum over scales of 20-30 A. Despite the eleven species present in the 
synthetic spectrum, it does not reproduce this small-scale structure. Branch et al. (2004) 
noted the existence of many weak features in the pre-maximum spectra of SN 2002cx. They 
were also unable to fit these features with Synow and speculated that either non-LTE level 
populations enhanced otherwise weak lines or clumpy velocity structure in the ejecta might 
be responsible. The SN 2005hk spectrum shows many of the same features, plus others 
visible due to its higher signal-to-noise ratio, making clumps at SN 2002cx-specific velocities 
seem less likely. In addition, if the ejecta were strongly clumped, we might expect to see 
some polarimetric signature, as in SN 2004dt (Wang et al. 2006; see below). We do not at 
this time have solid identifications for these features, but their confinement to wavelengths 
blueward of 6000 A suggests an origin in iron-peak elements. It is conceivable that many 
of them are also present in other SNe la, but are blended together in objects with higher 
expansion velocities. 



3.2. Spectropolarimetry 

3.2.1. Correction for ISP 

Before analyzing the polarimetry, we must correct the observations for the effects of 
interstellar polarization (ISP), both in our own Galaxy and in UGC 272. The expected 
Galactic contribution to the ISP is low, as the Galactic latitude is high (—63.5 deg) and the 
dust maps of Schlegel, Finkbeiner, & Davis (1998) predict a Galactic reddening, E(B — V), 
of only 0.023 mag. The polarization efficiency (pmax / E(B — V )) of Galactic dust has 
been found to be typically 3% mag -1 , with a maximum value of 9% mag -1 (Serkowski et al. 
1975). We observed two stars, HD 2158 and HD 2636, to act as probes of the actual Galactic 
ISP. Following the prescription of Tran (1995), both stars are within 0.6° of SN 2005hk and 
have spectroscopic parallaxes that place them more than 150 pc above the Galactic plane 
and thus sample almost all of the Galactic dust column along that line of sight. According 
to the discussion above, we expect a polarization of about 0.07%, with a maximum of 0.2%. 
The actual debiased, flux-weighted l^-band polarizations we measure for HD 2158 and HD 
2636 are 0.14% at 9 = 131° ± 2° and 0.25% at 9 = 130° ± 1°, respectively. 

Lacking a reason to prefer one star over the other, we averaged the two together to 
form our estimate of the Galactic contribution to the ISP, and smoothed versions of q and 
u were subtracted from the observations of the supernova. This non-parametric procedure 



- 9- 



was followed because at this low polarization level we must allow for the potential effects 
of instrumental polarization which might cause the polarization to not follow the standard 
wavelength dependence of Serkowski et al. (1975). The null standards, as noted above, had 
measured polarizations below 0.08%, but we have seen instrumental polarization artifacts 
at the 0.1% level with LRISp in the past, particularly on the blue side of the spectrograph 
(Leonard et al. 2002). The probe stars were observed immediately after SN 2005hk and with 
the slit at the same position angle, and thus should have similar amounts of instrumental 
polarization, if any. Any uncertainties in the removal of Galactic ISP are folded into the 
estimation of the host-galaxy contribution to ISP below. 

It is instructive at this point to examine the location of the SN 2005hk spectropolari- 
metric data points in the q — u plane, as shown in Figure 4. The data have been binned up 
to 100 observed A per point for clarity, and noisy data at the ends of the spectrum are not 
plotted in this figure. The first important fact is that the overall polarization is low, less 
than 0.5%. The second is that the points do not lie randomly in the plane. The distribution 
of points is elongated in a direction which does not intersect or point to the origin. This 
clearly means that host ISP cannot be the sole source of polarization in SN 2005hk (Howell 
et al. 2001; Leonard et al. 2002). The effect of vector subtraction of a small value of ISP 
is similar to changing the origin of the q — u plane, as the wavelength dependence of ISP 
is not strong (Serkowski et al. 1975). Furthermore, the data points at shorter wavelengths 
are systematically located below those at longer ones, indicating a gradient in polarization 
percentage in this object. 

The contribution of dust in the host galaxy of SN 2005hk, UGC 272, to the total ISP is 
harder to determine than that of our Galaxy. The magnitude of the expected polarization 
can be estimated as above from the reddening, and the reddening in turn can be estimated 
from the interstellar Na I D AA5890, 5896 absorption lines. Absorption due to Na I D 
can be seen in Figure 3 at the redshift of UGC 272, with an equivalent width (EW) of 
~0.35 A. Barbon et al. (1990) found a crude relationship between the reddening and the 
total equivalent width of Na I D absorption, E(B-V) = 0.25 mag A" 1 x EW(Na I D). This, 
combined with the Serkowski et al. (1975) relationship for polarization efficiency, yields an 
estimate for the typical expected host ISP of 0.26%. 

Munari & Zwitter (1997) have also derived a relationship for the reddening as a function 
of the EW of the Dl component of Na I absorption, but the two components of the host 
Na I D absorption are not resolved in our spectrum. The sightlines in the Munari & Zwitter 
(1997) study showed a variation in the ratio of EWs for Na I D1:D2 of 1.1:1 to 2:1. Using 
these values as a bracket for the true value, the range of expected polarizations is 0.18- 
0.24% for the Munari & Zwitter (1997) relationship, consistent with that predicted using 
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the Barbon et al. (1990) relationship. The dotted circle shown in Figure 4 has a radius of 
0.26% to show the magnitude of the expected host ISP relative to the data, though dust 
with maximal polarization efficiency could show three times as much polarization. 

The Na I D methods described above can only estimate the magnitude of host ISP, 
but not the direction. Fortunately, we can use the SN 2005hk polarimetry itself to estimate 
the host ISP. The optical spectrum of an SN la at blue wavelengths is expected to have 
almost no continuum polarization on theoretical grounds. Opacity due to many overlapping 
spectral features from iron-peak elements should act to depolarize the continuum, even if 
the continuum is intrinsically polarized by electron scattering in an aspherical atmosphere 
(Howell et al. 2001). Different authors in the literature have chosen different methods to take 
advantage of this assumption. Howell et al. (2001) used an ISP that made the polarization at 
the (arbitrary) blue end of their spectra almost zero. Wang et al. (2006) and Leonard et al. 
(2005) used the polarization at specific wavelengths of blue (A < 5000 A) emission features to 
define an ISP under the assumption that the observed polarization of these emission features 
is entirely due to ISP. We adopt a spectral window of 4000-4200 A in the SN rest frame 
to measure our host ISP. The exact endpoints of this window are arbitrary, but this region 
was chosen to avoid the strong Ca II and Fe III absorption lines at nearby wavelengths 
which could potentially show polarization features. Our measured value for (qisp,uisp) °f 
(0.009%, —0.272%) has a formal uncertainty of 0.018% in each Stokes parameter. This point 
is marked by a square in Figure 4, and is consistent with our expectations from the reddening 
arguments above. We have adopted a Serkowski et al. (1975) law with R v = 3.1 to subtract 
this ISP on the assumption that dust in UGC 272 resembles Galactic dust. 

As an aside, we note that it is possible to use these polarization measurements as 
a constraint on the reddening to SN 2005hk. One possible explanation for the low ap- 
parent luminosity of SN 2002cx-like objects is extinction, though this is hard to reconcile 
with the relatively blue early-time B — V color of SN 2002cx (Li et al. 2003). We have 
adopted a host ISP of 0.272%, which for standard dust polarization efficiencies corresponds 
to E(B — V) =0.091 mag. If this is combined with the Galactic contribution and Ry of 3.1 
is adopted, then we estimate a total extinction of Ay pa 0.35 mag, not nearly enough to 
conclude SN 2005hk was an SN la of intrinsically normal luminosity (Phillips et al. 2006, in 
prep.). 



3.2.2. SN 200 5hk Polarization 



The continuum polarization of SNe la can be understood as the result of electron scat- 
tering in an aspherical supernova atmosphere (Hoflich 1991; Jeffery 1991). If the projection 
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of the supernova on the sky were circularly symmetric, the local polarization vectors would 
cancel each other when averaged over the unresolved object and no polarization would be 
measured. Therefore, the continuum polarization level represents a meaningful measurement 
of the asphericity of the supernova. A number of authors have used the spectral region near 
7000 A to define a continuum polarization because of the relative lack of strong lines (Howell 
et al. 2001; Kasen et al. 2004; Leonard et al. 2005). Here we choose the (rest-frame) wave- 
length range 7000-7500 A to define a continuum region in SN 2005hk. The wavelengths were 
chosen to avoid the O I A7774 absorption trough. We measure Stokes parameters of (0.124%, 
0.088%) for this region prior to ISP subtraction, and mark this point with a diamond in Fig- 
ure 4. After subtraction of our chosen ISP, this continuum region has a polarization angle 
of 35.6° at a level of 0.36%, an entirely typical value for a normal SN la, as discussed below. 

The one-dimensional presentation of spectropolarimetric data is problematic because of 
the inherent two-dimensional nature of the measured Stokes parameters. In addition, at low 
polarizations the formal polarization (P = (q 2 + u 2 )^) is statistically biased. One solution 
is to rotate the q — u coordinate system by a convenient amount to new, rotated Stokes 
parameters (RSP) chosen to align one Stokes parameter (called qnsp here) with the axis 
of the system. In this new coordinate system, qnsp is an estimator of the total amount of 
polarization aligned with the axis of symmetry and ursp can show deviations from a single 
axis of symmetry in the system. 

For SN 2005hk, we have chosen to rotate the Stokes parameters by 35.6° to align q^sp 
with our chosen red continuum window of 7000-7500 A. These new rotated Stokes parameters 
are presented in Figure 5 and show a continuum polarization rising from zero near 4000 A 
(by construction for this choice of ISP) to ~0.4% at 9000 A. The flatness of the 6 curve and 
the lack of significant features in ursp show that a single axis of symmetry dominates the 
continuum polarization. Note that at very low polarization levels, such as seen at wavelengths 
less than 5000 A in our data, the angle of polarization becomes ill-defined. No strong 
line features are seen in the polarimetry, but there does appear to be a weak polarization 
modulation at the Fe III A5129 line. A 0.2% depression in qnsp indicates some depolarization 
and a smaller feature in ursp m ay indicate a rotation of the axis of symmetry. Other, smaller 
polarimetric features may also be present, but we caution against over- interpretation of weak 
(<0.1 %) spectropolarimetric features. 

4. Discussion 

Early attempts to measure broad-band polarization of SNe la were unsuccessful, both 
due to the low intrinsic polarization levels and the difficulty of disentangling the effects of 
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ISP, so intrinsic polarization remained undetected until SN 1996X showed marginal evidence 
at the ~0.3% level (Wang et al. 1997). Since then, intrinsic polarization has been detected 
in a number of SNe la, both in the continuum and in line features (Howell et al. 2001; Kasen 
et al. 2003; Wang et al. 2003; Leonard et al. 2005; Wang et al. 2006). 

There are no published spectropolarimetric data for SN 1991T-like or SN 2002cx-like 
objects to act as a basis for comparison with SN 2005hk, but there are a few relevant 
measurements. SN 2004bv had an early-time spectrum that showed prominent Fe III lines 
like SN 1991T (Foley et al. 2004), and Pereyra & Magalhaes (2005) obtained a single epoch 
of i?-band imaging polarimetry two weeks after maximum light. They found a polarization 
of 0.11% ± 0.02%, consistent with the polarization of the foreground stars in their field 
and therefore with low intrinsic supernova polarization. Leonard et al. (2000) reported on 
maximum-light spectropolarimetry of the unique SN la 2000cx, which resembled the SN 
1991T-like class of objects (Li et al. 2001b). They found a continuum polarization of ~0.5%, 
with a 0.3% modulation across the Si II line, but most of the continuum polarization may be 
due to ISP (D. Leonard, private com.). SNe la with normal spectra and luminosities have 
shown continuum polarizations at the level of ~0.3% (Wang et al. 1997, 2003; Leonard et al. 
2005). Two subluminous SNe la, SN 1997dt and SN 1999by (which resembled SN 1991bg), 
have shown somewhat higher continuum polarizations, up to 0.5% and 0.8%, respectively, 
with the exact values dependent on assumptions about ISP and the limited wavelength range 
of the observations for these two objects (Howell et al. 2001; Leonard et al. 2005). 

In summary, the magnitude of the continuum polarization in SN 2005hk is entirely con- 
sistent with observations of other SNe la. Core-collapse supernovae, by contrast, show po- 
larizations of up to a few percent (e.g., Leonard & Filippenko 2005). Current evidence seems 
to indicate that the core-collapse supernova explosion mechanism is fundamentally aspher- 
ical, and objects with lower polarization levels are those with extensive hydrogen envelopes 
that act to damp and hide asymmetries. The low polarization of thermonuclear supernovae 
indicates a more nearly spherical explosion mechanism. Whatever unknown mechanisms are 
generating asphericities in the explosions of normal SNe la, such as instabilities in the growth 
of a deflagration front (e.g., Ghezzi et al. 2004), are probably also responsible in SN 2005hk. 

The continuum polarization percentage in some other SNe la varies with wavelength. 
The most striking example is the low-luminosity SN 1999by (Howell et al. 2001), whose 
polarization rises dramatically to the red, as it does to a lesser extent in SN 2005hk. Thomson 
scattering is independent of wavelength, so the slope in polarization percentage is interpreted 
as a signature of the depolarizing effect of line opacity, which is stronger at UV and blue 
wavelengths than in the red (Howell et al. 2001). This result is dependent on the choice of ISP, 
which has been driven by the expectation of low intrinsic polarization at blue wavelengths. 
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SN 2001el provides a note of caution, as the ISP derived from the polarization of blue emission 
features at early times (Kasen et al. 2003) does not agree with the polarization at late times 
(Wang et al. 2003) when the intrinsic polarization should be low. We do not have multiple 
epochs of polarimetric data on SN 2005hk, so we cannot use the late-time polarization to 
test our choice of ISP. 

A number of SNe la have shown definite spectropolarimetric line features, but in 
SN 2005hk they are weak. The strongest line features (P <2.5%) have been seen in those 
SNe la that show high-velocity absorption at early times in the Ca IR triplet and the Si II 
lines. Kasen et al. (2003) showed that partial obscuration of the photosphere by clumps 
of intermediate-mass elements at high velocity provides an adequate description of the po- 
larimetric data for SN 200 lei, perhaps due to a gravitationally confined detonation (Kasen 
& Plewa 2004). Wang et al. (2006) preferred to explain the strong polarization features of 
SN 2004dt by clumpy protrusions of silicon into an oxygen-rich layer in the outer SN atmo- 
sphere, maybe as a result of turbulent burning processes during the explosion. SN 2005hk 
does not show high-velocity absorptions and the Si II and Ca II lines are weak, so these 
strong polarization modulations are not expected. We can at least conclude that the lack of 
strong polarization increases in the absorption troughs of SN 2005hk limits the presence of 
inhomogeneous clumps above the photosphere. 

SNe la lacking the high-velocity material show much smaller line polarization features, 
usually at the same level as the continuum polarization or less (Howell et al. 2001; Leonard 
et al. 2005). Howell et al. (2001) showed that, in general, strong absorption lines in SNe la 
should be associated with polarization minima as most photons have been absorbed and 
reemitted in the line, destroying the directional information imprinted by electron scatter- 
ing. The spectral features in SN 2005hk may be just too weak to leave strong polarization 
signatures. For example, some of the strongest polarization features have been seen in other 
SNe la in the Si II A6355 absorption line, but in SN 2005hk the absorption depth is only 
~12% of the continuum. In the simple aspherical toy models of Leonard & Filippenko (2001), 
such a feature would only be expected to result in a polarization feature that was at most 
~12% of the continuum polarization level. This model does not take into account optical 
depth effects, clumping, or multiple overlapping lines, but it may allow us to understand the 
lack of strong spectropolarimetric line features in SN 2005hk. For a continuum polarization 
of 0.4%, this model would predict a polarization modulation in the Si II line that is unde- 
tectable at our current signal-to-noise ratio. The strongest optical lines in SN 2005hk are 
the blue Fe III lines and we do see evidence for depolarization by 0.2% and maybe an angle 
rotation in A5129. The polarization signatures in the blue may be expected to be small (in 
the absence of clumping) if all the weak overlapping iron-peak lines destroy polarization of 
the pseudo-continuum at blue wavelengths, so the additional opacity of strong lines has little 



-14- 



additional depolarizing effect. 

Although no pre-existing supernova model explains all the peculiarities of SN 2002cx, 
one new model that does address some of the issues and makes specific polarimetric predic- 
tions for both the continuum and lines is the ejecta hole model of Kasen et al. (2004). They 
followed up on the simulations of Marietta et al. (2000), who determined that the donor star 
to the exploding white dwarf of an SN la could leave a noticeable imprint on the ejecta in 
the form of a hole. The opening angle of the hole was 30-40° for main sequence, sub-giant, 
or red giant donor stars. Kasen et al. (2004) found that sightlines down the hole allow an ob- 
server to see more deeply into the ejecta. This has several effects on the observed spectrum, 
including increasing the prominence of higher-ionization species, weakening absorption lines, 
and lowering the observed photospheric velocity relative to more typical sightlines. These 
effects seem promising for SNe 2002cx and 2005hk. Unfortunately, they also found that 
objects viewed down a hole would be somewhat brighter than normal. Kasen et al. avoided 
this difficulty by proposing that if SN 1991T-like objects are intrinsically normal SNe la 
viewed down a hole, then SN 2002cx-like objects could be the analogues of the subluminous 
SN 1991bg-like class of objects, but viewed down an ejecta hole. One difficulty with this 
scenario is that the SN 2002cx-like objects have only been seen in late-type spiral galaxies 
(Jha et al. 2006), while the SN 1991bg-like objects are more common in early-type galaxies 
(Howell 2001). 

With that in mind, Kasen et al. (2004) made specific spectropolarimetric predictions 
for the consequences of an ejecta hole. If the system were oriented such that our line of sight 
were aimed straight down the hole, then the projection onto the sky would be symmetrical 
and very little polarization should be expected. Such an exact alignment would be quite rare 
and even small angular offsets from the axis of symmetry should produce measurable effects. 
They made a sample calculation of a polarization spectrum for a line of sight oriented 
20° to the axis of their model hole (opening half- angle of 40°) and found that such an 
orientation gave very low continuum polarization (~0.1%, reflecting the fact that their model 
was symmetric apart from the hole) with strong line polarization peaks (~1% level). We 
measure a somewhat higher continuum polarization and no strong line features, in apparent 
contradiction to the expectations of this model. 



5. Conclusions 

We have presented single-epoch spectropolarimetry of SN 2005hk and an analysis of the 
high signal-to-noise-ratio total-flux spectrum. The flux spectrum is very similar to that of 
the peculiar SN la 2002cx. Both supernovae resemble the class of SN 1991T-like objects 
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near maximum light in that they show prominent lines of Fe III and only weak lines of 
intermediate-mass elements like Si II and Ca II that characterize the early spectra of normal 
SNe la (Filippenko 1997). In addition, SN 2005hk shows low expansion velocities of ~6000 
km s -1 , similar to those of SN 2002cx. Li et al. (2003) found that no existing theoretical 
model for SNe la could reproduce all of the peculiarities of SN 2002cx. Observations of 
SN 2005hk provide additional constraints for models of this new subclass of SNe la. 

We have presented a set of line identifications using Synow. Our Synow fit for 
SN 2005hk is very similar to a previous analysis of SN 2002cx (Branch et al. 2004), but 
with five additional ionic species. Four of these are due to our superior wavelength range 
and do not necessarily indicate spectral differences between the two objects. A detailed 
analysis of a time series of spectra could also clarify the identification of spectral features 
attributed to those species whose presence is less certain in our analysis of this individual 
spectrum. If all of our possible line identifications are correct, then we have identified the 
signatures of iron-peak elements, intermediate-mass elements, and maybe unburned carbon 
at similar velocities in the early spectra of SN 2005hk. Our SYNOW-based analysis is un- 
suitable for determining the actual relative abundances of these three groups of species, but 
a clear next step for future research is for detailed supernova spectrum synthesis codes to 
attempt to model the spectrum of SN 2005hk. 

Of particular interest is whether normal deflagration products (e.g., Nomoto et al. 1984) 
can fit the SN 2005hk spectrum for favorable values of the temperature, density, and ion- 
ization. Jha et al. (2006) found evidence for the large-scale mixing predicted in three- 
dimensional simulations of deflagrations (Reinecke, Hillebrandt, & Niemeyer 2002; Gamezo 
et al. 2003) in the presence of low-velocity intermediate-mass elements and possibly O I at 
late times in SN 2002cx. Alternatively, our identification of Co II and Ni II in SN 2005hk 
might be evidence for freshly synthesized 56 Ni and its decay products in the outer layers of 
the supernova, as was interpreted in SN 1991T to be the sign of a detonation (Filippenko et 
al. 1992a; Mazzali et al. 1995). The short half-life of 56 Ni implies that nickel must have been 
present in larger quantities at earlier times for it to still be seen in the current spectrum, 
so the signatures of nickel should be looked for in any spectra of SN 2002cx-like objects 
obtained at very early epochs. A potential consequence is that radioactive material present 
in the outermost ejecta might not make a large contribution to the luminosity. 

One of the most perplexing mysteries of SN 2002cx is how that object could have a 
spectrum dominated by iron at both early and late times and yet have such a low luminosity 
and small expansion velocities. The characteristic velocity scale (and hence kinetic energy 
per unit mass) of SNe la near maximum light (>10,000 km s _1 ) is set by the energy release 
of the fusion of carbon and oxygen up to iron-peak or even just intermediate-mass elements. 
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The presence of a large fraction of unburned material that adds inertial mass to the ejecta 
could solve this puzzle, but evidence of this raw material in the spectra is scarce. One 
argument against a large amount of unburned material is the rapid spectral evolution seen 
originally in SN 2002cx, and confirmed in SN 2005hk by the stronger lines of intermediate- 
mass elements seen in this spectrum compared to the one from three days earlier (Jha et al. 
2006). This may indicate that the outer ejecta have relatively low mass, though an analysis 
of the photometry and spectral evolution is needed to confirm this hypothesis. 

One line of evidence in favor of unburned material is the possible presence of weak O I 
lines in the late-time spectra of SN 2002cx. Here we have looked for carbon in the early- 
time spectrum of SN 2005hk and have obtained some potential evidence for C III. Features 
suggestive of C II are also present, but we were unable to find a good fit, though more 
detailed handling of non-LTE effects might allow a secure identification to be made. Similar 
features seen in normal SNe la have been identified as C II in the past (Branch et al. 2003; 
Garavini et al. 2004). Hoflich et al. (2002) were able to show that the peculiar subluminous 
SN 1999by had unburned material above the photosphere at early times by detecting C I 
lines in the near-IR. Near-IR spectroscopy of SN 2005hk or another member of the class 
of SN 2002cx-like objects at early times would be very useful to see if they more closely 
resemble SN 1999by in this respect or the more normal SN la sample of Marion et al. (2006), 
which lacked the C I lines. 

Our spectropolarimetric data for SN 2005hk are perhaps the least unusual aspect of the 
object. They show a low level of continuum polarization (~0.4%) and evidence for only weak 
line features. Therefore, the level of asphericity is quite normal for an SN la and models for 
the unusual aspects of the SN 2002cx-like class of objects that invoke large deviations from 
spherical symmetry can be rejected, such as the ejecta hole model of Kasen et al. (2004). 
The lack of strong line features in the spectropolarimetry may be related to the weakness 
of the Si II and Ca II lines. If this is correct, SN 1991T-like objects observed in the future 
may only show weak polarimetric line features, unless clumping in the ejecta is important. 
Lastly, a full identification of the weak spectral features might provide additional insight into 
the nature of SN 2005hk. 
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Table 1. Synow Fitting Parameters 



C III 1 Mg II Si II Si III S II Ca II Ti II Fe III Co II Ni II 

X ref (A) 4647 7773 4481 6347 4553 5454 3934 4550 4420 4161 4067 
t 0.3 0.08 0.35 0.2 0.3 0.25 0.8 0.04 0.7 0.2 0.1 
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Fig. 1. — Flux spectrum of SN 2005hk compared to the similar objects SN 2002cx (2002 May 
17; Li et al. 2003), SN 1997br (1997 April 16; Li et al. 1999), and SN 1991T (1991 April 
19; Filippenko et al. 1992a). Also plotted is a representative normal SN la, SN 1994D (1994 
March 17; Filippenko 1997). Narrow nebular emission lines from a superposed H II region 
have been clipped from the SN 2002cx spectrum. Numbers in parentheses are dates relative 
to S-band maximum. The strong Fe III absorption lines are clearly more blueshifted in SN 
1991T and SN 1997br than in SN 2005hk. Some strong absorption lines due to intermediate- 
mass elements (S II, Si II, and Ca II) are marked on the SN 1994D spectrum, but they are 
weak or absent in the other four objects. For complete line identifications, see Figure 2. 
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Fig. 2. — Synow fits (dashed lines) to the flux spectrum of SN 2005hk (solid line). The 
fit presented in the top panel has seven ions we believe are present with high confidence. 
The fit shown in the bottom panel is the same as that in the top, but with four additional 
species whose presence is less certain (Co II, Ni II, Ti II, and C III). The flux scale is linear 
and approximates f v with an arbitrary tilt applied to better show spectral features and the 
wavelength scale is logarithmic to better show relative velocity widths of spectral features. 
The primary species responsible for spectral features are marked. See the text and Table 1 
for Synow fit parameters. 
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Fig. 3. — Section of the SN 2005hk spectrum enlarged to better show the numerous uniden- 
tified weak spectral features. The weak interstellar Na I D lines from our Galaxy (MW) and 
the host galaxy of SN 2005hk are labeled. 



-25 - 




-0.4 -0.2 0.0 0.2 0.4 0.6 

q(%) 

Fig. 4. — Polarization data for SN 2005hk (after correction for Galactic ISP) are plotted in 
the q—u plane. The data have been binned to 100 A for clarity, and only points with observed 
wavelengths between 4000 A and 9000 A are shown here. The symbol size is proportional to 
wavelength, so bins at longer wavelengths are represented by larger points. Bluer-wavelength 
bins tend to be found in the cloud of points at lower values of u than redder ones. The 
dotted circle represents the typical magnitude of ISP expected from the host galaxy based 
on analysis of the Na I D lines. The square marks the location of our chosen ISP vector, 
and the diamond marks the location of our continuum measurement (7000-7500 A). The 
line between them represents the direction of the new q axis after subtracting the ISP vector 
and rotating the Stokes parameters about the new origin at the square. See text for details 
and Figure 5 for the spectropolarimetry after ISP subtraction and rotation. The cross at 
(0.45,-0.35) shows the median statistical error bars for the points in the plot. 
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Fig. 5. — Polarization data for SN 2005hk after removal of ISP and rotation about the new 
origin to align the rotated Stokes parameter qRSP with the continuum axis of symmetry. 
The flux spectrum (scale is linear fx) is overplotted on the qnsp data. The polarization data 
are presented in 50 A bins. A continuum polarization of 0.4% can be seen in addition to a 
feature at the Fe III line near 5000 A. Any polarization features in ursp would represent 
deviations from a single axis of symmetry. The relatively flat 9 curve over those wavelengths 
where the polarization angle is numerically well-defined confirms that SN 2005hk has a 
single dominant axis of symmetry. The small gap in the data near 5600 A is the region of 
the spectrum affected by the dichroic and where small spurious polarization features were 
present. 



